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NOMENCLATURE 


Symbol 

Dimensionless  Form 

Definition 

at 

Effective  aspect  ratio  of  tail  appendag 

A 

At 

- a/i2 

Area  of  tail. 

I , I 

y Vo 

Moment  of  inertia  of  body  and  of  dis- 
placed fluid  about  the  y-axis. 

kr 

k 2’  k' 

Hydrodynamic  mass  coefficients. 

i 

V 

= 1 

Characteristic  length  of  submarine. 

m 

o 

m ' 
o 

= mo4pi3 

Mass  of  submarine. 

m2 

= mo/(i+k2)/IPr 

Mass  coefficient. 

M,  M2 

M ' 

= M/-^pi3V2 

Hydrodynamic  moments  about  the 
y-axis  through  the  c.g. 

M 

q 

M 

w 

Md 

M 1 

q 

M ' 
w 

MI 

0 

* m,4m4v 
- mw4"<3v  •{ 

= m54p<3v2 

Derivatives  of  moment  component 
with  respect  to  angular  velocity 
component  q,  velocity  component  w, 
and  stern  plane  angle  6 . 

n ’ 

y 

= (I  +1  k')/ipi5 

Moment  of  inertia  coefficient  about 
the  y-axis 

q 

q' 

= x2(s)  = qjf/V 

Angular  velocity  component  relative 
to  y-axis. 

t 

s 

= tv/i 

T ime 

u,  V,  w, 

Velocity  components  of  origin  of 
body  axes  relative  to  fluid. 

V 

Velocity  of  origin  of  body  relative 
to  the  fluid. 

Z f Z ^ ) Z 2 

Z' 

= z/^pi  V 

Hydrodynamic  normal  force,  positive 
downward. 
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Symbol 


L 

q 


z 

w 


Z6 


NOMENCLATURE 

(continued) 


Dimensionless  Form 


Z ' = Z /ip|3V 

q q 2 

Zw  3 Zw4P,2v 

z6  = Z6  4p'ZV2 


< 


Definition 


Derivatives  of  force  component  with 
respect  to  angular  velocity  com- 
ponent q,  velocity  component  w, 
and  stern  plane  angle  5 . 


^t.ail 


See  Equations  5 and  6,  and  text. 


z 


o 


Depth  of  c.g.  from  horizontal  ref- 
erence line. 


a 


a = Xj(s) 


The  angle  of  attack. 


a 


tail 


The  local  angle  of  attack  at  the  tail. 


«3S  = a a 


Angular  displacement  of  stern 
planes,  positive  trailing  edge  down. 


9 


9 


Angle  of  pitch  or  inclination  the 
body  axis  from  the  horizontal. 


P 


P* 


Mass  density  of  water. 
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INTRODUCTION 

The  theoretical  problem  of  predicting  the  course  which  a ship  follows 
in  response  to  a prescribed  stern-plane  or  rudder  motion  has  not  been 
satisfactorily  solved,  A satisfactory  solution  of  this  problem  would  have 
considerable  practical  value,  particularly  for  the  motion  of  a submarine 
where  "mane uve rability  in  depth**  is  so  important.  The  theory  is  unsatis- 
factory in  the  sense  that  trajectories  describing  the  motion  cannot  reliably 
be  predicted  from  the  results  of  captive  model  tests.  The  commonly  ac- 
cepted reason  for  this  failure  is  that  the  hydrodynamic  characteristics  which 
are  determined  from  captive  model  tests  are  not  known  with  sufficient  pre- 
cision or  reliability.  This  reason  is  quite  plausible  in  view  of  the  discrep- 
ant results  which  are  often  obtained  from  different  model  tests  of  the  same 
prototype . 

The  present  study  was  undertaken  with  the  object  of  ascertaining 
whether  it  is  possible  to  combine  information  from  captive  model  tests 
with  information  from  free -running  model  tests  for  the  purpose  ~f  constuct- 
mg  hydrodynamic  characteristics  which  may  be  introduced  into  differential 
equations  of  a given  form,  such  equations  being  used  to  characterize  the 
motion,  and  thus  allowing  one  to  predict  new  trajectories.  Only  hand  compu- 
tation methods  have  been  used  in  this  investigation. 

Although  no  very  fixed  conclusions  can  be  drawn,  it  appears  that  more 
accurate  free-running  tests  (as  well  as  captive  model  tests)  are  required  fo 
obtain  positive  results.  However,  it  is  misleading  to  imply  that  nothing  of 
value  can  be  learncid  from  this  type  of  approach.  One  can,  in  fact,  roughly 
predict  some  trajectories  and,  given  enough  patience,  one  could,  perhaps, 
continue  to  modify  the  equations  so  that  they  fit  more  and  more  trajectories. 
It  is  known  that  nonlinear  differential  equations  must  be  used  to  character- 
ize the  motion  and  it  is,  therefore,  evident  a priori  that  theoretically  there 
must  always  exist  some  ambiguity  regarding  the  validity  of  the  equations. 
However,  from  a practical  point  of  view  one  may  say  that  if  equations  have 
been  constructed  which  have  as  one  solution  a given  trajectory,  then  these 
equations  should  be  approximately  valid  for  trajectories  which  are  not  too 
different,  i.e.  , trajectories  which  do  not  involve  a different  type  of  man- 
euver and  do  not  involve  large  differences  in  the  magnitudes  of  any  of  the 
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parameters  which  directly  influence  the  motion. 

GENERAL  THEORY 

A great  deal  has  been  written  on  the  theory  of  the  motion  of  a solid 
body  through  a fluid  and  here,  as  elsewhere,  the  usual  assumptions  of  mo- 
tion in  undisturbed  water  ard  in  a plane  (in  this  case  vertical)  are  made. 

The  coordinates  describing  the  motion,  of  the  center  of  gravity  are  taken 
along,  and  perpendicular  to,  the  longitudinal  axis  of  the  ship.  This  is  cus- 
tomarily done  because  the  hydrodynamic  inertia  forces,  which  depend  on 
the  geometry  of  the  body,  are  most  simply  expressed  in  this  form.  We  are 
interested  in  a description  of  the  motion  in  a direction  perpendicular  to  the 
main  line  of  motion  as  the  ship,  being  self-propelled,  does  not  accelerate 
or  decelerate  appreciably  along  the  trajectory  while  turning.  The  viscous 
forces,  which  affect  flow  mainly  in  the  boundary  layer,  are  therefore  neg- 
lected as  they  have  little  influence  on  the  motion  except  insofar  as  they 
determine  ihe  propeller  force  required  to  maintain  a constant  speed.  The 
coordinate  system,  terminology,  and  notation  follows  that  prescribed  and 
set  forth  in  Ref.  1.  The  general  equations  of  motion  have  been  derived  in 
Refs.  2,  3,  and  4,  and  a complete  derivation  will  not  be  given  here.  How- 
ever, some  description  of  some  of  the  various  forms  the  equations  may  take 
is  necessary  to  understand  the  significance  of  the  empirically  determined 
hydrodynamic  forces. 

Equations  of  Motion 

The  classical  equations  of  motion  for  a body  of  revolution  with  inertia 
coefficients  k^,  and  k1,  moving  in  one  plane  through  an  ideal  fluid  are, 
(see  Ref.  2) 

m (1  + k, ) u + m (1+k-J  w H = X,, 

o 1 ov  2'  1 

m (1  + k^)  w - mQ(l  + k^)  u 0 = ZV  (1) 

I (1-fk)  9 + m (k,  * k,)  uw  = M, , 

y ' ' O 1 2 1 

where  X^,  Zj,  and  represent  external  and/or  hydrodynamic  (or  hydro- 
static) forces  arising  in  part  from  the  nonideal  character  of  the  fluid. 

There  may  be  included  in  Xj  , Zj,  and  forces  arising  from  the  rudder. 
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stern  plane,  or  the  tail  fins,  and  gravity  or  buoyancy  forces. 

If  the  body  is  self-propelled,  and  has  not  toe  large  an  angle  of  attack, 
then  the  nu"  component  of  the  velocity  is  considered  to  be  essentially  con- 
stant, and  the  first  equation  contains  only  second  order  terms  and  is  elimi- 
nated. This  assumption  results  in  two  equations  for  w and  G . 

It  may  be  noted  at  this  point  that  the  dependent  variable  w is  some- 
times replaced  by  the  angle  of  attack  a . When  the  latter  is  small  and 
when  the  speed  V is  essentially  constant,  the  relations  u =*  V,  w r Vq  , 
allow  a simple  transformation  of  the  dependent  variables. 


The  mass  and  the  inertia  1^  are  quite  unambiguous.  The  inertia 

coefficients  k^,  k^,  and  k1  may,  in  theory,  be  calculated  from  the  geometry 
of  the  body,  and  ordinarily  an  estimate  is  obtained  in  this  way.  However, 
the  hydrodynamic  forces  represented  by  these  coefficients  are  only  a part 
of  the  total  hydrodynamic  forces  and  sometimes  the  equations  are  written 
in  a form  such  that  on  the  left  hand  side  one  has  only  the  body  inertia  forces 
while  on  the  right  hand  side  one  has  the  total  hydrodynamic  (and  other)  forces, 

e.  g.  , 


m(w-uG)  = Z, 
ov 

I 0 = M 

y 


(2) 


The  hydrodynamic  forces  and  moments,  Z and  M,  will  depend  not  only  on 
w anc  0 bat  also  on  the  "accelerations"  w and  6,  Furthermore,  as  has 
already  been  pointed  out,  there  will  be  included  in  Z and  M the  stern  plane 
and  fin  effects  as  well  as  hydrostatic  forces,  consequently  some  discussion 
of  these  forces  is  in  order. 

Forces  and  Moments 

The  mathematical  formulation  of  the  problem  reduces  to  a detailed 
specification  of  what  is  to  be  included  in,  and  what  excluded  from,  the  terras 
Z and  M. 


In  the  first  place,  one  may  point  out  that  besides  the  five  variables,  w, 
0,  w,  9,  and  u,  there  is  the  stern  plane  angle  ^(t),  and  that,  strictly 
speaking,  Z and  M depend  functionally  on  all  of  these  quantities  in  a very 
complex  nonlinear  way.  Several  simplifications  of  the  problem  are  necessary 
and  the  Trst  of  these  is  the  assumption  that  only  linear  terms  in  the 
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accelerations  w and  0 are  involved  and  that  the  principal  ones,  of  these 

terms,  have  the  values  predicted,  e.  e.  . Z = ~k->m  , Mx  = - k‘I  , whereas 

r & w Z o 0 y 

the  coupling  terms  Zk  and  M.  are  negligible . A second  simplification  re- 

w v/  u ^ 

suits  from  the  assumption  that  all  hydrodynamic  forces  are  proportional  to  u , 
v/ith  the  consequence  that,  except  for  hydrostatic  forces  or  moments,  the  equa- 
tions may  be  made  independent  of  the  speed.  The  present  discussion  relates  to 
neutrally  buoyant  dynamically  stable  bodies  with  a small  metacentric  height 
and  traveling  at  high  speed  so  that  hydrostatic  forces  and  moments  are  ne- 
glected relative  to  hydrodynamic  ones. 

With  these  simplifications  the  functions  Z^  and  are  functions 

only  of  the  three  variables  w,  0,  and  6g  = 6,  the  first  two  depending  im- 
plicitly, the  latter  explicitly,  on  time.  The  subscript  s on  6 will  hence- 
forth be  omitted  as  this  should  cause  no  confusion. 


It  is  known  that  within  the  range  of  variables  which  occurs,  the  func- 
tions Zj  and  Mj  are  "nonlinear M in  all  three  variables.  However,  in 
view  of  the  explicit  dependence  of  6 on  t it  is  considered  that  the  prob- 
lem is  analytically  unapproac Viable  without  a linearization  with  respect  to 
6 . Such  a linearization  may,  of  course,  be  done  in  several  different  ways. 
In  the  present  paper  the  initial  slope  is  taken  as  the  basis  for  the  linear- 
ization, To  summarize  then,  we  have 


Z = Z,  - m k^  w + m k.  u 0 , 

I O L O 1 

M - M , -m  (k,  - k-> ) uw  - k1  I 0 , 

1 ox  1 2 y 

and, 

Z j = Z ^ (w,  0,  <5 ) = Zg  6 + (w  , 9)  , 
(w,  9,6)  = *Mg  6 + M^(w,  9)  . 


(3) 


(4) 


Z2  an<3  M2  are  nonlinear  functions  of  w and  0,  and  Zg  and  Mg  are  stern 
plane  coefficients  obtained  by  an  appropriate  linearization. 


HYDRODYNAMIC  CHARACTERISTICS 

With  the  assumptions  of  the  preceding  section,  it  is  seen  that  the  im- 
portant hydrodynamic  characteristics  which  must  be  determined  are 
Z Jw,  0)  and  M-,(w,  0),  as  well  as  the  stern  plane  coefficients  Zg  and  Mg  . 
The  latter  depend  on  how  the  equations  are  linearized  with  respect  to  6 . 
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The  lift  Z(6),  and  moment  M(6),  associated  with  a given  stern  plane 
angle  6 may  be  determined  for  zero  angle  of  attack,  or  transverse  velocity 
w,  and  for  zero  turning  rate  9.  When  so  determined,  the  slopes  dZ(6)/dfi 
and  ,dM(6)/dC  are  often  found  to  be  quite  constant  throughout  a large  range 
of  6 (Ref.  8)  and  the  values  of  these  slopes  at  6 = 0 are  the  stern-plane 

coefficients  Zg  and  Mg  . If,  however,  the  angle  of  attack  (for  example) 
is  not  zero,  and  similar  measurements  are  made,  it  is  found  that  the  slopes 
Zq  and  Mg  are  not  essentially  constant  throughout  the  region  of  interest. 
One  may,  as  has  been  done  in  Ref.  5,  average  the  coefficients  over  a speci- 
fied range  of  the  variables  w and  0.  It  is  less  arbitrary,  although  perhaps 
also  less  satisfactor  y,  to  determine  Zg  and  Mg  by  the  first  method  and 
this  is  what  has  been  done  in  the  present  treatment. 


The  characteristics  Z^(w,  0)  and  M^(w,  0)  are  nonlinear  functions  of 
\a  and  0 and  may  be  expressed  giaphically  or  as  polynomials  in  w and  0 . 
The  general  character  of  the  functions  seems  well  known,  although  there 
seems  to  be  some  doubt  about  the  reproducibility  of  the  results  from  one 
measurement  to  another.  Furthermore,  in  the  present  instance  there  are 
incomplete  experimental  data  available  relating  to  the  specific  tail  configur- 
ation of  interest.  A certain  amount  of  effort  has  been  directed  towards  the 
prediction  of  the  linear  terms  of  these  functions.,  e.g.,  Refs.  6 and  1 . It 
is  pointed  out  in  Ref.  6 that  a crude  but  fairly  reliable  simplification  of  the 
problem  results  from  the  assumption  that  the  lift  Zj  can  be  split  into  two 
components,  as  has  been  done  in  Eq.  4,  one  due  to  the  deflection  of  the 
stern,  plane  and  one  due  to  the  tail  fins  which  act  as  lifting  surfaces.  The 
latter  component  will  be  proportional  to  the  local  angle  o£  attack  at  the  tail 
and  can  be  estimated  using  a result  of  aerodynamic  lifting  theory.  The  first 
component  will  be  proportional  to  the  stern  plane  angle  and  must  be  deter- 
mined empirically.  Thus  one  may  write  (see,  tor  example,  Ref.  6), 


5Z,  + a t 

0 tail  tail 


w he  re , 


i nd , 


' tail 


2ffat\ 

+ 2 


piV 


Effective -aspect- ratio  of  tail, 


(5) 

(6) 


A - Dimensionless  area  of  tail. 
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Furthermore,  under  the  simplifying  assumption  that  the  dimensionless  dis- 
tance from  the  c.  g.  to  the  tail  is  one-half,  then  the  angle  of  attack  at  the 
tail  is  just, 

<*4,  = W/V  + M/2V.  (7) 


Thus, 


Zx  = 5 Z5  + wZt/v  + 9iZt/2V  , 


Similarly,  the  moment  is  taken  to  be, 

M,  = 6M£  + (wZt/2V  + 6iZt/4V)  f 


(8) 

(9) 


since  the  factor  l/ 2 represents  the  moment  arm. 


Except  then  for  the  two  empirically  determined  coefficients  Zg  and 
ivig  , one  could  in  principle  calculate  all  terras  occurring  in  the  two  linear- 
ized differential  equations  from  a knowledge  of  the  geometry  of  the  body 
alone.  In  practice,  however,  these  results  are  not  sufficiently  reliable  and 
use  is  made  of  information  obtained  from  hydrodynamic  dynamometer  tests. 
In  order  to  relate  this  information  to  the  above  estimated  coefficients,  it  is 
necessary  to  note  that  what  is  measured  in  the  captive  model  tests  includes 
all  hydrodynamic  terms  but  does  not  include  the  acceleration  terms  relating 
to  the  body  alone. 


When  Eqs.  (8)  and  (9)  are  introduced  into  Eqs.  (3)  and  (4),  the  follow- 
ing results  are  obtained: 

Z = - m k2w  + 6 Zg  + w2t/v  + ©(m^u  + IZJtV), 

(10) 

M * - k'ly  9 + 6Mfi  + wjzt</2V  - m^kj-  k2)  u]  + 9f2 Zt/4V  . 
Therefore, 


M 3 m u(k,-  k.)  + izjzv  , 

w o x Z 1 ' t 

= mokt  u + fZt/2V  , 

Mg  = f2 Zt/4V  . 


(U) 
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In  dimer  sionle  ss  notation  and  tailing  ~~  0,  zrr : k'  — 1 these  results 

take  the  form 


'"W  . I 


w 


2tt  a-tAt 

a + 2 


M ’ = l Z ' + m ’ , 
w 2 w o 


(12) 


Z ’ = iz1  , 

q 2 w 


M ’ = i z 1 
q 4 w 


The  foregoing  considerations  should  be  viewed  as  suggestive  rather 
than  determinative  and  if  tbe  free -running  tests  appear  to  be  more  consis- 
tent with  slightly  different  linearized  coefficients,  such  coefficients  would 
have  just  as  much  logical  justification  as  the  "predicted”  coefficients  of 
Eq.  (12). 

As  a starting  point,  however,  it  is  convenient  to  calculate  coefficients 
using  Eq.  (12)  to  modify  these  in  a consistent  way  on  the  basis  of  captive 
model  tests,  and  then  perhaps  to  modify  them  again  on  the  basis  of  free- 
running  te  sts . 

The  application  of  these  considerations  may  now  be  made  to  the  tests 
of  the  free -running  model  of  the  Hydrodynamics  Laboratory.  The  model  is 
scaled  down  1:100  from  the  SST,  Scheme  IV  submarine.  No  hydrodynamic 
tests  have  been  made  on  the  model  and  all  hydrodynamic  tests  which  have 
been  made  using  captive  models  of  the  prototype  are  based  on  tail  assemblies 
that  differ  to  some  extent  from  the  tail  of  the  Hydrodynamics  Laboratory 
model. 

The  alternate  tail  "B"  (for  which  test  results  are  available)  and  the 
model  tail  may  be  considered  to  be  roughly  similar  in  size  and  shape  and 
one  might  expect  that,  although  the  model  tail  has  a somewhat  larger  pro- 
jected area,  the  experimentally  determined  coefficients  for  the  alternate 
tail  ,!B"  could  be  used  to  estimate  the  coefficients  tor  the  model  tail.  The 
question  arises  as  to  what  sort  of  correction  factor  should  be  applied. 

Insofar  as  it  does  not  se^m  to  be  possible  to  precisely  calculate  coefficients 
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from  a knowledge  cf  the  geometry,  the  method  employed  results  in  an  esti- 
mate only. 

The  result  of  applying  the  simplified  theory  to  the  alternate  tail  "B'\ 
using  a projected  area  of  299*4  ft  and  an  aspect  ratio  of  2.  28,  is  shown  in 
Table  1,  which  also  contains  the  tabulated  experimental  results  on  tail  "B". 
The  approximate  ratio  of  the  calculated  results  to  the  experimental  results 
for  tail  "B"  is  taken  as  a correction  factor  by  means  of  which  ’’predicted" 
coefficients  for  the  model  tail  are  obtained  from  calculated  coefficients.  A 
third  set  of  coefficients  labeled  "revised"  is  shown  and  will  be  explained 
later.  Stern  plane  coefficients  are  also  shown  in  Table  1.  For  tail  "B" 
these  are  the  measured  coefficients.  For  the  model  tail  these  are  calcu- 
lated from  the  steady-state  turning  results  by  the  method  outlined  in  a later 
section  of  the  paper. 


TABLE  1. 


Coefficients  for  the  SST,  Scheme  IV. 


Alternate  Tail  "B" 

Correction 

Factor 

Hydrodynamics  Lab.  Model  Tail 

Calculated 

Measured* 

Calculated 

Predicted 

Revised 

z ' 

w 

- . 0249 

- . 0250 

1. 00 

-.  0277 

-.0277 

-.  02832 

4.  Vi 

w 

j_  nncoc 

T • v w ^ ; w 

+ . 0075 

n on 
v . u v 

■f.  0044 

+. 0055 

. A a / *»  r 

+ . UUOL? 

Zq 

- . 01245 

- . 0083 

1 . 50 

01385 

-. 00923 

-. 00863 

Mq 

- . 00622 

- .0065 

1.00 

-.  00693 

00693 

-.00625 

Z6 

- 

- . 0070 

- 

- 

-.  0090 

-. 00697 

M6 

- 

! - . 0040 

- 

“ 

-.  00514 

-.00412 

The  measured  coefiicients  on  tail  "B"  were  obtained  by  estimating  the 
slopes  (at  the  origin)  of  the  hydrodynamic  characteristics  given  in  Ref.  7 . 
The  values  of  and  Mg  are  contained  in  Ref.  8. 


The  preceding  discussion  related  to  the  linear  terms  of  the  nonlinear 
functions  Z^(v/,  9)  and  M^(w,  9).  The  nonlinear  terms  have  in  some  cases 
been  determined  experimentally  and  application  of  these  terms  to  the  dif- 
ferential equations  has  been  made.  e.  g.  , Ref.  5.  However,  for  the  model 
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configuration  these  nonlinearities  can  only  be  inferred  from  a knowledge  of 
the  nonlinearities  for  a closely  similar  configuration.  In  a later  section  of 
this  paper  a method  of  constructing  these  nonlinearities  from  the  trajectories 
of  the  free -running  model  is  indicated. 

The  stability  coefficients  (at  infinite  speed),  <Tj  and  have  been  cal- 

culated using  the  "revised"  coefficients.  They  have  the  values  «j\  = - 0.465, 
"z  = -3.835. 


ANALYSIS 

Solution  of  the  Linearized  Equations 

The  preceding  section  has  shown  how  all  terms  may  be  estimated  for 
the  linearized  equations.  These  were  not,  however,  put  in  dimensionless 
form.  When  written  in  the  customary  dimensionless  form  they  are: 

nij  da/ds  - u.  Z - q ' (m  1 + Z 1 ) ~ Z 1 6 (s)  , 

L a o q O 

(13) 

n ’ dq/ds  - a.  M ' - q’  M ' = M ^ 6 ( s)  , 
y ^ a q o 

where  a and  6 are  measured  from  the  values  of  a and  6 corresponding  to 
a steady,  horizontal,  neutrally  buoyant  run.  All  the  coefficients  Z^,  etc.  , 
are  measured  for  this  value  of  a and  for  q'  = 0.  For  convenience,  the 
equations  have  been  rewritten  in  the  form, 


dxj /ds 

+ + a12x2  = cx  (s), 

dx^/ds 

+ a21Xj  + a,2x2  = c2  (s), 

(14) 

where  x^(s) 

= a (s),  and  x2(s)  ~ q'  (s). 

These  have  been  solved  and,  also,  9(s)  = 
6(s)  of  the  form, 

r 

j q'(s)ds  has  been  computed  for 

6(s)  = 

J 6os/ T » 0 — * — t . 

| <5 

[_  o > s > T • 

(15) 

The  solutions  involve  the  stability  exponents  (at  infinite  speed)  cr^  and  , 
as  well  as  various  other  constants  derived  from  the  a’s  and  c's.  The  solu- 
tion is  straightforward  and  only  the  results  corresponding  to  the  initial 
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conditions,  x^(0)  = x2(0)  = 0 will  be  given  here 


(6, 


h C,  (e  1 - <r.  s - 1)  + D,  (« 

L 2 2 


2J  . 


a,  s - 1) 


, 0 ^ S ^ T, 


X.  (s)  = < 
*2  ' 


X1  (T)  + (6n /T) 

*2  ° 


<T,  T O'.(S-T) 

C , ( e - 1)  ( e - l)  + D,  (e 

L ll  l2 


V 


<r?(s  - x) 

(e  - - 1) 


1) 

(16) 


(6q/t)  j^C2A^(c  1 - o-j  s - 1)  + {l>z/*z)  ( 


+ IB2s2 


(s)  = ^ 


J 


s > T , 


a2s 

e - s - 1 ) 


, 0 ^ s ^ T, 


(17) 


e7(a)  + B26o(s  - t)  + (djx) 


<r.T  «r,(s  - x) 

(C2/«r1)(e  1 - 1 )( e 1 - 1) 


(TyX  a (s  - t) 

C * \ / £ 


+ (D2/a2)(e  - 1 ) ( e 


1) 


, s y x . 


The  stability  exponents  are3 


°1,2  = “ 2 ^all  + a22^  -"\/  (all  + a22^  ' * ” iall  a22  ’ a12a21  ^ * (18) 


The  constants  which  enter  into  the  solution  ar; 


C1  *^a22  C1  ”a12  c£  ^ /°i 
*2  _H  1 2 JjA  1 

(»i  - ^ 


D, 


c 1 "^(a?2  C1  ”a12  C 2 ^ ^2 
L2  2 ^21  1 


(19) 


^2  <°'2  * V 


B, 


( a22  C1  * a12  c2  ^ ^^l  a 2 
^11  2 ^21  4 
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If  the  initial  conditions  £x^(0)  and  x^O)”]  , are  not  zero,  correction  terms 
must  be  added;  however,  in  the  model  tests  any  deviation  from  zero  initial 
conditions  is  experimentally  not  detectable  and  all  the  calculations  assume 
Xj(0)  = X£(0)  = 0.  The  steady-s4ate  values  of  Xj(s)  and  x^s)  *or  a con" 
stant  stern-plane  angle  5 are  easily  shown  to  be 


x,  (oo)  = B 
1 


l2  ° • 


(20) 


The  sveady-state  turning  rate  x^oo)  = q'(oo)  is  related  to  the  steady-state 
turning  radius  R(oo)  by  virtue  of  the  relation 


I ! 


1 / \ i /-»  / \ 

q (oo)  = i/n\oo;  , 


(21) 


where  R is  measured  in  ship  lengths 


* ! 
l 

cyi 


Treatment  of  the  Nonlinear  Equations 

If  the  functions  and.  are  not  linearized,  then  the  dimensionless 

equations  corf e sponding  to  Eqs.  (|3)  might  be  written  in  the  form 


da/d 


s - 


Z'z(a,  q*)  4-  m^q'J/ rn^  = Zq  6(s)/m . 


(22) 


dq'/ds  - q'  )/n  ' = Mfi'  6 (s)/n 


In  the  nonlinear  form  of  th^se  equations,  the  functions  J^Z  ^(a,  q')  + m^q^j/m^ 
and  M^(a,  q ) may  be  represented  fjraphically  as  functions  of  a for  different 
values  of  q*  . 

In  the  vicinity  of  a = q'  - c functions  will  be  straight  lines 

whose  slopes  and  intercepts  are  determined  by  the  coefficients  Z^, 
and  . The  trajectories  of  the  model  are  used  to  estimate  the  shape  of 
the  hydrodynamic  characteristics  when  a and  q’  are  not  small.  It  is  1 ike  Iy 
that  trajectories  which  arc  predicted  from  characteristics  constructed  in 
this  way  will  only  be  approximate  and  that  the  approximation  is  less  likely  to 
be  satisfactory  the  more  the  predicted  trajectory  differs  from  those  used  in 
constructing  the  hydrodynamic  characteristics. 

APPLICATION  TO  THE  MODEL 

Several  sets  of  calculations  have  been  carried  out  to  construct  the 

CONFIDENTS  L 


*T’ 

r* 


CONFIDENTIAL 


-12- 


hvdr  odynarn.ic  characteristics  and  M2  ir>  such  a way  that  trajectories 
calculated  from  them  will  give  results  which  are  consistent  with  the  experi- 
mental results. 

It  is  characteristic  of  the  free -flight  tests  that  there  is  no  very  sensi- 
tive test  of  the  validity  of  the  linear  terms  of  these  functions.  This  is,  per- 
haps, mainly  due  to  the  technical  limitations  of  the  free -flight  tests,  i.  e.  , 
trajectories  which  one  might  conceive  of  as  being  always  within  the  linear 
range,  and  preferably  steady-state  trajectories,  are  not  possible  within 
the  limited  space  of  the  testing  tank.  On  the  other  hand,  those  trajectories 
which  are  possible,  and  which  have  reasonably  small  scatter,  have  only  a 
very  limited  "linear"  range  which  is  not  adequate  to  test  the  predicted  or 
assumed  linear  terms. 

It  is  also  characteristic  of  the  free -flight  test*  lhat  the  two  dependent 
variables  of  the  differential  equations,  namely  the  angle  attack  a and  the 
angular  turning  rate  q\  are  not  measured  directly  but  can  only  be  derived 
from  ofher  measured  quantities.  The  turning  rate  q'  is  the  derivative  of 
the  body  inclination  0 which  one  can  measure,  and  the  angle  of  attack  is 
even  less  precisely  known  since  it  is  measured  by  the  difference  between 
the  body  inclination  0 , and  the  inclination  of  the  trajectory  tan  * dz/ds. 

The  scatter  in  the  values  of  a and  q1  from  presumably  identical  trajectories 
will  therefore  be  much  greater  than  the  scatter  in  the  values  for  the  inclina- 
tion 0 and  the  depth  z. 

It  might  justifiably  be  argued  that  one  is  only  interested  in  the  in- 
clination and  the  depth,  and  that  if  differential  equations  can  be  constructed 
which  enable  one  to  predict  these  quantities,  th  n one  need  not  be  concerned 
about  the  logical  basis  for  the  equations.  One  must,  however,  always  keep 
in  mind  the  danger  of  extrapolation  from  such  results.  It  has  been  pointed 
out  that  the  "predicted"  coefficients  given  in  Table  1 for  the  model  are  not 
predicated  by  any  very  compelling  reasons.  On  the  other  hand,  the  data 
obtained  from  the  trajectories  are  insufficient  to  allow  one  to  derive  these 
coefficients.  As  a consequence,  the  selection  of  an  appropriate  set  of 
linear  terms  for  the  characteristics  is  reduced  to  a trial  and  error  proce- 
dure. Calculations  (of  a type  to  be  described)  were  made  using  the  "pre- 
dicted" coefficients  and  using  several  sets  of  revised  coefficients.  The 
"revised"  coefficients  shown  in  Table  1 are  therefore  somewhat  arbitrary. 
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However,  they  lead  to  results  which,  within  a limited  range,  appear  to  be 
quite  reasonable.  Only  the  results  of  calculations  using  the  "revised"  coef- 
ficients will  be  presented  here. 

Figure  1 shows  the  steady-state  turning  rate,  q1  (oo)  and  angle  of  at- 
tack, a(oo)  for  a constant  stern  plane  angle  6^  = <5^.  The  dashed  curve  is 
the  experimental  curve  obtained  with  the  model  tor  dives.  The  straight 
lines  correspond  to  Eq.  (20).  Numerically, 

q'(oo)  = - 1.42  , and 

a (oo)  = **  0.  820  <5 
' ' o 


The  first  value  is  obtained  by  approximating  the  slope  of  the  experimental 
curve.  The  latter  value  is  obtained  by  taking,  somewhat  arbitrarily,  the 
value  0.577  for  the  ratio  Z^1  /Mg.  The  solid  curves  of  Fig.  1 are  smooth 
curves  constructed  to  approach  the  straight  lines  in  the  linear  region.  In 
the  case  of  the  curve  for  q'(oo)  the  curve  is  constructed  to  approximate  the 
experimental  curve.  In  the  case  of  the  curve  for  a(oo)  the  experimental 
data  are  not  shown  because  of  the  large  scatter  and  the  curve  shown  is  some- 
what arbitrary  because  of  this  scatter. 


The  procedure  whereby  the  characteristics  j^Z^a,  q' ) + q'J  /m-J 
and  M^(a,  q')/n^  shown  in  Figs.  ? and  3 are  constructed  will  now  be  outlined. 


Using  the  "revised1*  coefficients  given  in  Table  1,  together  with  the 
values  = 0.01840,  ml,  = 0.0354  and  n^  = 0.001787,  the  linear  portion  of 
the  characteristics  becomes, 


(Zl,  + m ? q'  )/m\  = - 0.  800  a + 0.  290  qr  , and 

^ D 6 


M!,/V  = + 3.  50q  - 3. 50  q'  . 

Furthermore,  zl/m'^  = - 0.1969  and  Mg/ny  = “2.  307.  Thus,  Eqs.  (22)  be- 
come 


da/ds  + 0.  800a  - 0.  290  q*  + nonlinear  terms  = - 0.  1969  6(s), 
d4/ds  - 3.  50  a + 3.  50  q'  + nonlinear  terms  - - 2.  307  6(s). 


(23) 


For  steady-state  turning  dq'/ds  = da/ds  = 0,  and  for  each  value  of  6 (s)  = 6q 
there  are  values  of  q'(oo)  ar\d  a(oo)  taken  from  Fig.  1 which,  through  Eqs.  (23), 
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determine  the  steady-state  curves  in  Figs.  2 and  3.  The  appropriate  values 
of  q*  are  indicated  on  the  steady-state  curve. 

The  dive  maneuver  resulting  from  a stern  plane  deflection  of  18°, 
achieved  at  a 5°  per  second  rate,  ha9  been  used  as  a second  set  of  data  from 
which  Z2  and  Mp  may  be  constructed.  Figure  4,  which  shows  the  inclina- 
tion 0(s),  and  Fig.  5,  which  shows  the  depth  zq(s)  include  experimental  and 
calculated  curves  for  incline  runs  as  well  as  the  dive.  However,  only  one 
of  these,  e.  g.  , SST^  Run  S-101  has  been  used  in  constructing  the  charac- 
teristics. Figure  4 also  shows  the  stern  plane  action  for  the  dive  and  in- 
cline trajectories.  The  short  vertical  lines  which  intersect  the  trajectories 
in  Figs.  2 through  7 may  be  correlated  with  the  stern  plane  action. 

The  desired  information  which  is  to  be  obtained  from  these  curves 
consists  of  the  angular  velocity  q'(s)  = d9/ds,  the  angular  acceleration 
dq'/ds,  the  angle  of  attack  a (s)  and  the  rate  of  change  of  a,  da/ds.  This 
information,  which  is  shown  in  Figs.  6 and  7,  is  most  satisfactorily  ob- 
tained by  graphical  differentiation  with  subsequent  numerical  integration  as 
a check.  Other  methods,  including  difference  tables  and  the  fitting  of  high 
degree  polynomials  to  0(s)  and  zq(s)  with  subsequent  differentiation  have 
been  tried  but  do  not  seem  to  the  writer  to  be  any  more  satisfactory. 

It  is  found  in  practice,  as  explained  elsewhere,  that  the  initial  por- 
tion of  the  curves  cannot  be  obtained  with  satisfactory  accuracy  by  this 
procedure.  Consequently,  Eqs.  (16)  and  (17)  with  the  appropriate  numeri- 
cal constants  obtained  from  the  "revised”  column  of  Table  1,  have  been 
used  to  compute  the  initial  portions  of  the  "18°  dive  trajectories"  in  Figs. 

2 and  3.  The  trajectories  are  continued  using  the  values  of  da/ds,  dq'/ds, 
a,  and  q'  from  Figs.  6 and  7 in  Eqs.  (22).  As  before,  the  values  of  q'  are 
marked  on  the  trajectory. 

A first  approximation  to  the  hydrodynamic  characteristics  may  be 
obtained  by  drawing  smooth  curves  through  corresponding  valuer,  of  q'  on 
the  steady-state  curve  and  the  18°  dive  curve.  Logically,  it  would  seem 
to  be  more  reasonable  to  repeat  this  procedure  for  the  other  dives  (15°, 

12°,  9U»  etc.  ) and  in  this  way  to  obtain  several  points  on  the  characteristic 
curve  corresponding  to  each  value  of  a'.  However,  the  data  are  not  quite 
accurate  enough  to  make  this  procedure  as  satisfactory  as  it  seems  to  be 
to  construct  the  smooth  curves  by  judgment  based  on  similar  experimentally 
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obtained  hydrodynamic  characteristics. 

Having  constructed  the  hydrodynamic  characteristics,  it  is  then  pos- 
sible to  compare  a new  computed  trajectory  with  one  obtained  experimentally. 
For  this  comparison,  the  incline  trajectory  obtained  by  throwing  the  stern- 
plane  to  18°  at  5°/sec,  holding  for  4.  8 seconds  (prototype),  and  returning  to 
neutral  at  5°  sec,  was  used.  The  computed  trajectory  is  shown,  together 
with  two  of  the  experimental  trajectories  SST^  Runs  196  and  197,  in  Figs. 

4 and  5.  The  computed  trajectory  was  obtained  by  a step-by-step  numerical 
integration  using  values  of  and  obtained  from  Figs.  Z and  3.  Experi- 
mentally it  is  difficult  to  run  an  incline  trajectory  which  exactly  repeats  the 
initial  portion  of  the  dive  trajectory  as  it  theoretically  should  do;  consequently, 
one  cannot  expect  to  obtain  much  better  agreement  than  that  shown. 

The  logical  extension  of  this  procedure  would  be  to  calculate  a "pullout" 
trajectory  using  the  incline  trajectory  up  to  s » Z.  56,  and  to  continue  the 
calculations  with  a reverse  stern  plane  action  for  a prescribed  time  duration. 
One  would,  of  course,  have  to  extrapolate  the  hydrodynamic  characteristics 
in  some  way. 

This  has  been  attempted;  however  the  method  has  a practical  limita- 
tion since  the  numerically  integrated  solution  tends  to  diverge  from  the  true 
solution  and  if  this  divergence  becomes  too  large  the  "solution"  begins  to 
oscillate  and  no  longer  gives  a true  prediction.  This  could  perhaps  be 
avoided  at  some  expense  of  effort;  however,  because  of  the  arbitrariness 
involved  in  the  extrapolation  of  the  hydrodynamic  characteristics,  it  does 
not  seem  to  be  justified. 
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SUMMARY  AND  CONCLUSIONS 

The  paper  represents  an  attempt  to  make  a theoretical  analysis  of 
free -running  trajectories  of  a model  experimental  submarine,  the  SST, 
Scheme  IV.  The  basic  assumptions  on  which  the  theory  is  based  are  set 
forth  in  some  detail.  The  equations  of  motion  are  not  derived,  this  aspect 
of  the  problem  having  been  adequately  covered  elsewhere.  However,  the 
various  forms  which  the  equations  may  take  are  indicated  and  explained, 
particularly  with  reference  to  the  hydrodynamic  forces  and  moments  in- 
volved. These  hydrodynamic  terms  are  discussed  separately  at  some 
length  to  form  a basis  for  the  estimation  of  numerical  values  for  the  model. 
Finally,  the  analysis  consists  essentially  in  the  construction  of  hydrodynamic 
characteristics  using  available  information  from  both  static  tests  and  free- 
running  tests.  The  use  of  these  characteristics  to  predict  new  trajectories 
is  demonstrated. 

There  is  at  present  rather  too  much  scatter  in  the  data  from  the  free- 
running  tests  to  allow  of  its  being  systematically  used  in  the  way  in  which  it 
has  been  used  in  the  analysis  reported  here.  However,  with  more  accurate 
data  it  does  seem  possible  that  more  extended  results  of  the  form  reported 
here  could  be  obtained  and  that  such  data  might  have  value  for  design  pur- 
pose s . 
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